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Abstract 

By molecular dynamic simulations we study a system of particles interacting through a contin- 
uous isotropic pairwise core-softened potential consisting of a repulsive shoulder and an attractive 
well. The model displays a phase diagram with three fluid phases, a gas-liquid critical point, a 
liquid-liquid critical point, and anomalies in density, diffusion and structure. The hierarchy of the 
anomalies is the same as for water. Here we study in a systematic way the effect on the anomalies 
of varying the softness of the potential. We find that, making the soft-core steeper and more 
penetrable, the regions of density and diffusion anomalies contract in the T — p plane, while the 
region of structural anomaly is weakly affected. Therefore, a liquid can have anomalous structural 
behavior without having density or diffusion anomalies. We show that, by considering as effective 
distances those corresponding to the maxima of the first two peaks of the radial distribution func- 
tion g(r) in the high-density liquid, we can generalize to continuous two-scales potentials a criterion 
for the occurrence of the anomalies of density and diffusion, originally proposed for discontinuous 
potentials. However, we observe that the knowledge of the structural behavior within the first 
two coordination shells of the liquid is not enough to establish, in general, the occurrence of the 
anomalies. By introducing the density derivative of the the cumulative order integral of the excess 
entropy, measuring shell by shell the amount of order in the liquid, we show that the anomalous 
behavior is regulated by the structural order at distances as large as the fourth coordination shell. 
By comparing the results for different softness of the potential, we conclude that the disappear- 
ing of the density and diffusion anomalies for the steeper potentials is due to a more structured 
short-range order. All these results increase our understanding on how, knowing the interaction 
potential, we can evaluate the possible presence of anomalies for a liquid. 

PACS numbers: 64.70.Ja, 82.70.Dd, 61.20.Ja, 64.70.qj, 65.20.De 
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I. INTRODUCTION 



Several liquids exhibit anomalous behaviours, i.e. behaviours that differ from what is 
expected for normal liquids such as argon. Water is the most common example of such 
systems. For example, while normal liquids contract when they are cooled, water expands 
under T = 4°C at ambient pressure^. 

However, water is not unique in this respect. The anomlaous behavior in density, in fact, 
has been observed also in experimets for Bi£, Gar-, Te^, Ge^Te^ and simulations for 
silica^— , silicon^- and BeF^-. All these systems, including water, present a temperature 
of maximum density (TMD) below which density decreases when temperature is lowered at 
constant pressure. 

Nevertheless, liquid water is anomalous also in other properties; for example, in its dy- 
namics and its structure. While the diffuson coefficient D of a normal liquid decreases 
as density or pressure are increased, water is characterized by a region of the phase dia- 
gram where D increases when the pressure is increased. Experiments show that the normal 
behaviour is restored at pressures higher than P ~ 1.1 kbar at 283K -. 

Regarding the structure, normal liquids tend to become more structurated when com- 
pressed. This can be quantified by two order parameters: a translational order parameter 
t, quantifying the tendency of pairs of molecules to adopt preferential separations, and an 
orientational order parameter Qi that measures the tendency of a molecule and its nearest 
neigb hours to assume a specific local arrangement, as considered by Steinhardt et al.— . For 
a normal liquid t and Qi increase with pressure. Water, instead, shows a region where the 
structural order parameters decrease for increasing pressure or density at constant T, i.e. 
the system becomes more disordered, as shown by Errington and Debendetti with molecular 
dynamics simualtionsiS. Also simualtions for a model of silica display similar anomalies, as 
shown by Shell et al.-. 

The sequence of the anomalies in water is well determined in experiments 1 and 
simulations^. For the T — p phase diagram the water structural anomaly region is en- 
compassing the diffusion anomaly region which includes the density anomaly region. 

Despite the fact that these anomalies can be found also in other liquids the sequence 
of anomalies may be different or may be the same. It is different, for example, for silica 
where the anomalous difussion region contains the structural anomaly regions that, in turn, 
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includes the density anomaly region^. The sequence of anomalies is, instead, the same for 
several isotropic models 14 . 

Isotropic models are systems of identical particles interacting through a central potential. 
It has been shown that isotropic potentials with a repulsive soft-core, either a ramp or a 
shoulder, have anomalies that resemble those of water—"—. 

Here we consider one of these potentials showing anomalies: the Continous Shoudered 
Well potential (CSW), introduced in^ 4 - and studied in detail in^ 4 -. This model can be con- 
sidered as the continuous version of the Discontinuous Shouldered Well potential (DSW) 
studied in 4 ^— . The two potentials have been shown to share the same qualitative phase 
diagram, with a liquid-iquid critical point in the supercooled liquid phase. However, while 
the CSW has anomalies, the DSW has not. The aim of the present study is to understand 
the origin of this difference. In particular, we tune the CSW potential in such a way to 
approximate, at least around the repulsive shoulder, the DSW, changing in this way the 
softness of the potential. Our results show that the continuous tuning of the CSW induces a 
continuous change on the regions of density and diffusion anomaly, but not a relevant effect 
on the structural anomaly region. When the CSW approches the limit of the DSW, the 
regions of density and diffusion anomaly shrink as they would disappear. 

The paper is organized as follows. We present the model in Sec. 2 and give details about 
the simulations in Sec 3. In Sec 4 we present our results about the phase diagram, and in 
Sec 5 our results about the anomalies, first with a criteria based on the excess entropy and 
then with a direct calculation. In Sec 6 we discuss the results, we generalize some proposed 
criteria and we study the effect of the long-range order. In Sec 7 we give our conclusions. 



II. THE MODEL 

We study the CSW model^ 4 - consisting of a set of identical particles interacting through 
the isotropic pairwise potential given by 



U(r) 



U, 



R 



— Ua exp 



25\ 



+ 



21 



(1) 



1 + exp (A(r-R R )/a 

where a is the diameter of the particles, Ra and Rr are the distance of the attractive 
minimum and the repulsive radius, respectivelly, Ua and Ur are the energies of the attractive 
well and the repulsive shoulder, respectivelly, 8\ is the variance of the Gaussian centered in 
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FIG. 1: Continuous lines represent the CSW potential for different values of the parameter A (slope 
at r = Rr). Dotted black line represents the discountinuous shouldered well potential (DSW)^~— . 
By increasing the values of A the CSW potential tends to DSW around Rr. 

Ra, and A is the parameter which controls the slope between the shoulder and the well at 
-Rr (Fig. [1]). Varying the parameters the potential can be tuned from a repulsive shoulder 
to a deep double well. In particular, by increasing A the soft-core repulsion becomes more 
penetrable near the minimum of the attractive well, and the softness of the potential increases 
for r > R R and decreases for r < R R . 

The set of values chosen for our analysis are as in^^, Ur/Ua = 2, Rr/ci = 1.6, Ra/ci = 2, 
(5a/ a) 2 = 0.1, while we tune A. In order to reduce the computational effort, we establish 
a cuttoff for U(r) at r c /a = 3, because U(r) — >■ for r/a > 3 for all the considered values 
of A. In particular, we consider the values A = 15, 30, 100, 300, 500 (Fig. [1]) going from the 
case A = 15 studied in Ref.— ,— to slopes that aproach the infinite value of DSW. 

For A = 15 this model is known to present anomalies in density, diffusion and 
structure^ 1 ^. Its phase diagram displays a liquid-gas phase transition and a liquid-liquid 
phase transition, both ending in critical points^ 1 ^. 
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III. MOLECULAR DYNAMICS 



For every chosen value of A we simulate the behaviour of the system governed by the 
CSW potential using standard molecular dynamics techniques 4 ^. We perform simulations 
in the NVT ensemble for N = 1372 particles interacting in a cubic box of volume V with 
periodic boundary conditions at temperature T. We use the Allen thermostat keeping T 
constant by reescaling the velocities of the particles at each time step by a factor (T/T)*, 
where T is the instantaneous kinetic temperature and T the fixed temperature of the thermal 
bath 4 -. We integrate the equations of motion using velocity Verlet integrator—. By fixing p 
and T, we calculate the pressure in terms of the second virial coefficient 46 . All calculations 
are presented in dimensionless reduced units in terms of the particle diameter and the depth 
of the attractive well: P* = Pa 3 /U A , T* = k B T/U A , p* = pa 3 , and D* = D{m/a 2 U A ) 1/2 , 
where D is the diffusion coefficient. 

We use a simulation time step 5t* = 3.2 x 10~ 3 defined in units of (a 2 m/U a) 1 ^ 2 (of the 
order of ~ 1.7 x 10~ 12 s for water-like molecules and « 2.1 x 10~ 12 s for argon-like atoms). 
Simulations are equilibrated after 10 4 time steps, followed by 10 6 time steps during which 
thermodynamic quantities are computed every 100 time steps. Positions and velocities are 
stored every 1000 time steps for the calculations of structure and dynamic quantities. For 
each value of p* and T* we average over two runs, starting from independent configurations. 

Near the liquid-liquid critical point (described in the next section) we observe that the 
liquid phase is metastable with respect to the crystal phase. The nucleation of the crystal 
is marked by a subsequent drop in the potential energey 43 -. The lifetime of the metastable 
phase is defined as the time between the end of the equilibration and the energy drop. 

IV. THE PHASE DIAGRAM 

First we verify that our simulations reproduce the known phase diagram for the model 
with A = 15^^. Next, we calculate the phase diagram for the other considered values of 
A, finding always the same qualitative behaviour (Figj2]). In particular, we find two regions 
with coexisting fluid phases: (i) gas and liquid at low p*, (ii) low density liquid (LDL) and 
high density liquid (HDL) at higher p*. The two coexistence regions end in critical points: 
C\ for gas-liquid coexistence, C2 for LDL-HDL coexistence (Table HJ. For A < 100 the HDL 
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FIG. 2: Isotherms in the P* — p* phase diagram from simulations for CSW potential with A = 30 
(from top to bottom, T* = 1.4, 1.3, 1.2, 1.1, 1.0, 0.9, 0.8, 0.7, 0.6, 0.55, 0.5 and 0.45). At low 
T*, isotherms have no monothonic behaviour and show two van der Waals loops, corresponing to 
a liquid-gas (low densities) and a liquid-liquid (higher densities) first order phase transitions both 
ending in critical points (C\ and C2 respectivelly). Inset: the crossing of the isotherms signals the 
presence of the density anomally. In this and the following figures errors, where not shown, are 
smaller than symbol size. 

phase is metastable with respect to the crystal, but with a lifetime long enough to allow us 
to equilibrate the liquid around C2. For A > 100 the HDL lifetime is too short to equilibrate 
the liquid around C2 and we extrapolate the location of C2 from higher T* isotherms. 

As A is increased both C\ and C2 tend to the corresponding values for the DSW potential. 
This result is consistent with the idea that the DSW can be seen as a limiting case of the 
familly of CSW potentials studied here. However, temperature and pressure of the LDL- 
HDL critical point for A = 500 are still far from the values for the DSW potential. 

V. HIERARCHY OF THE ANOMALIES FOR DIFFERENT A 

Barros de Oliveira et al.— show that the CSW potential has anomalies in density, diffusion 
and structure. Furthermore, they found a hierarchy of the anomalies as reported for other 
two-scale potentials and for the SPC/E water. Here we study in a systematic way how these 
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p* 


pK 


T 6 2 


p* 

r c 2 


P*C 2 


A = 


15 


0.95 ±0.06 


0.019 ±0.008 


0.08 ± 0.03 


0.49 ±0.01 


0.285 ±0.007 


0.247 ± 0.008 


A = 


30 


1.01 ±0.07 


0.022 ± 0.008 


0.08 ± 0.03 


0.50 ±0.01 


0.267 ± 0.005 


0.255 ± 0.009 


A = 


100 


1.06 ±0.04 


0.025 ± 0.005 


0.08 ± 0.03 


0.53 ±0.02 


0.243 ± 0.005 


0.262 ± 0.008 


A = 


300 


1.06 ±0.05 


0.027 ± 0.009 


0.09 ± 0.02 


0.51 ±0.01 


0.231 ±0.007 


0.271 ± 0.008 


A = 


500 


1.08 ±0.06 


0.027 ± 0.008 


0.09 ± 0.03 


0.52 ±0.01 


0.204 ±0.007 


0.272 ± 0.008 


DSW 


1.24 ±0.01 


0.030 ±0.001 


0.09 ± 0.02 


0.69 ± 0.02 


0.110 ±0.002 


0.280 ± 0.020 



TABLE I: Temperatures Tq and Tq , pressures P c and Pq , densities p c and p c (in reduced 
units) corresponding to the gas-liquid critical points C\ and the LDL-HDL critical point C2, respec- 
tively, for different values of A for the CSW potential. We estimate the critical parameters by using 
the condition (d 2 P/dp 2 ) \t = at C\ and C2. We confirm the data for A = 15 from Ref.— «. For 
sake of comparison, we indicate also the values for the DSW potential with parameters Ur/Ua = 2, 
R R /a = 1.6, R A /a = 2, and (5 A /a) 2 = 0.1 from Ref.i^ 

anomalies depend on the parameter A of the CSW potential. In Fig. [3] we summarize all 
our results. In the next subsections we discuss the details of Fig. |3j 

A. Excess entropy 

As a starting point for the exploration of the anomalous regions of the potential and in 
order to predict whether a given anomaly will appear, we perform an analysis of the excess 
entropy. The excess entropy s ex is denned as the difference between the entropy s of a real 
fluid at a given T and p and the entropy of an ideal gas s lg under the same conditions, 
s ex = s — s lg , and quantifies the entropy lost due to interactions and correlations between 
the particles of the fluid. It has been shown that the excess entropy allows an approximate 
estimate of the regions where density, diffusion and structure are anomalous. In particular, 
Errington et al.— proposed that the condition under which the density anomaly appears 
is given by E ex > 1, where by definition is S cx = (ds ex / d\a p)x- The diffusion anomaly 
can also be predicted, following Rosenfeld's empirical parametrization^, by the condition 
S ex > 0.42. Structural anomaly is determined by the criterion S ex > 0, because for normal 
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FIG. 3: Hierarchy of anomalies (as denned in Sec. V) in the T* — p* plane, plotted for (a) A = 15, 
(b) A = 30, (c) A = 100 and (d) A = 500. As the value of A increases, both the regions of 
density anomaly (•) and diffusion anomaly (■) contract with the diffusion anomaly region always 
encompassing the TMD line. The region of structural anomaly (between A-symbols for Qq maxima 
and the high-p* T-symbols for t minima) is only weakly affected. 

fluids the excess entropy is always decreasing when density increases isothermically^. 

The evaluation of the excess entropy migth be computationally expensive. However, it 
can be aproximated directly by the two-body contribution term 



which depends only on the radial distribution function and represents the dominant contri- 
bution to the excess entropy^— (between 85% and 95% of s ex for a Lennard- Jones fluid^ 1 ^). 
Since the integrand of Eq.(J2]) goes to zero for large r, i.e. where g{r) goes to 1, we calculate 
s 2 up to a cutoff given by half the simulation-box size, distance at which we check that 
g(r) fa 1 within the range of T and p considered in this work. From Eq.(J2]) we define and 
calculate E 2 = (<9s 2 /<91n We observe that the predicted regions of anomalies get narrow 
as the slope A increases (Fig. H]). 

Since s 2 < s ex , then it is always S cx < E 2 . Therefore, the criteria E 2 > 1, E 2 > 0.42 and 
E 2 > overestimate the regions of density, diffusion and structural anomalies, respectively, 
with about 30% uncertainty, according to Rosenfeld 48 . As we will discuss in Section VI, we 
observe that the anomalous regions predicted in Fig. H] are indeed an overstimate of the one 
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FIG. 4: Excess entropy S2 (upper panels) and its derivative £2 (lower panels) for the potentials 
corresponding to A = 30 (left panels) and A = 500 (right panels), for T* = 0.5 (O), 0.6 (□), 0.7 
(0), 0.8 (a), 0.9 (v). Horizontal solid lines in the bottom panels are the thresholds for £2 for the 
anomalies (from top to bottom for density, diffusion and structural anomaly, respectivelly) . As 
described in the text, these criteria overestimate the regions of anomalies. 

we find as a result of the detailed analysis presented in the next sections. Nevertheless, the 
approximate criteria based on E 2 give a preliminary idea of where the anomaly occur. 

B. The density anomaly 

For all the considered values of A, we find a temperature of minimum pressure along 
the isochores near the LDL-HDL critical point C 2 (Fig. EJ). These temperatures correspond 
to the TMD line at constant P. We find that the anomalous region where p decreases for 
decreasing T shrinks for increasing A and possibly tends to collapse onto one single point 
in the P* — T* plane for A — > 00. This result would be consistent with the behaviour of the 
DSW potential, that does not show density anomaly^"—. 
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Temperature T 

FIG. 5: The P* — T* phase diagram for A = 30 (a) and A = 500 (b). In both panels, we show 
isochores (with constant density separation 5p* = 0.01) from p* = 0.04 (at high T* and low P*) 
to 0.20, from 0.205 to 0.215 (with 5p* = 0.005), and (with dp* = 0.01) from 0.22 to 0.25 (at low 
T* and high P*). The points where isochores cross correspond to the gas-liquid and LDL-HDL 
critical points (C% and C2, respectivelly) . The black dashed line is a guide for the eye estimating 
the TMD line, correspondig to the line of minima along the isochores. 

C. Diffusion anomaly 

For the family of potentials given by Eq. (1) with the considered values of A we find 
an anomalous-diffusion region, i.e. a region (pomin < P < PDmax) where D increases with 
increasing density at constant T (Fig. H}. The diffusion coeficient is calculated from the 
mean square displacement of a single particle 

(Ar(t) 2 )^([r(t + t)-r(t )] 2 ), (3) 



11 




Density p* 



FIG. 6: Diffusion coefficient D* as a function of density p* for A = 30 and isotherms from T* = 0.45 
(bottom) to T* = 1.9 (top). Inset: Isotherms from T* = 0.45 (bottom) to T* = 0.7 (top) Black 
solid lines delimit the anomalous region where D* grows from -D m in to D max for increasing density 
between p* n < p* < p* n 

where t Q is any time at equilibrium and the average is over the initial t a and over the particles 
in the system, taking the long-time limit 

n .. <Ar(i) 2 > 

D = hm -V (4) 

t^oc 2dt v ' 

where d = 3 is the dimension of the system. We find that as in the case of the density 
anomaly, the region of diffusion anomaly contracts as A increases. We find that the diffusion 
anomaly region always encompasses the TMD line between pDmin and pDmax (Fig- E]). 



D. Structural anomaly 

The structural behaviour of the system is studied by calculating the translational or- 
der parameter t and the orientational order parameter Q$ defined in the following. The 
translational order parameter is defined as*^ 1 ^ 

POO 

t= / |<7(0-1|^ (5) 

where £ = rp 1 / 3 is a reduced distance (in units of the mean interparticle separation p -1 ^ 3 ) 
and is the radial distribution function. Since in our simulations g(£) ^ 1 for £ > £ c = 8, 
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FIG. 7: Translational order parameter t as a function of density p* for A = 30 and isotherms from 
T* = 0.5 (top) to T* = 1.9 (bottom). Dotted lines joining the maxima and minima along the 
isotherms are a guide for the eye delimiting the region between pt max and pt min , where t decreases 
for increasing p*. 

we integrate up to a cuttoff distance £ c . As the parameter t depends only on the deviations 
of <?(£) from unity, its value is sensible to long range periodicities. For an ideal gas g(£) 
is constant and equal to 1 and there is no translational order (t = 0). For a crystal phase 
g(£) 7^ 1 for long distances and t becomes large. While normal fluids show a monotonic 
increase of t with increasing density, for the CSW potential t increases for p < pt max , and 
reaches a maximum at pt max . Above Pt max , for increasing p, t decreases until it reaches a 
minimum at pt min . For p > pt min , t recovers the normal behaviour (Fig. [7]). 
The local orientational order for a particle i and an index / is defined as^ 



average of the spherical harmonics Yj m with indices I and m, evaluated over the vectorial 
distance between particles i and j. 

The calculations involved in the evaluation of Q\ can be computationally expensive. In 
order to reduce the computational effort, we express the spherical harmonics in terms of the 





where k is a fixed number of nearest neighbour particles and {Y^f, = 1 Ylj=i ^im( r «i) i s the 



13 




FIG. 8: Orientational order parameter Qq as a function of p* for A = 30, for isotherms going from 
T* = 0.5 (top) to T* = 1.9 (bottom). The dotted line is a guide for the eye joining the maxima 
along the isotherms. 



Legendre Polynomials Pi obtaining 

1/2 

(7) 



2 k-x 



3=1 h=j+l 



where 7* h is the angle formed by the vectors and r ih . 

A global order parameter is obtained by averaging the local order Q\ over all particles 

1 N 

G'=Af£<#- (8) 

i=l 

We study the case I = 6 and k = 12 as in Ref.— . For I = 6 the sixth order Legendre Poly- 
nomial is given by P%(x) = ^ (231x 6 — 315x 4 + 105x 2 — 5). For an isotropic homogeneous 
system, Q e reaches its minimum value Qf = 1/y/k = 0.287 (k = 12), while for a fully 
ordered f.c.c arrane ement Q% c = 0.574. 

The normal beahviour of Qq for a fluid without anomalies is to increase monotonically 
when p* increases at constant T*, or when T* decreases at constant p*. Qq decreases to the 
value Q'f at high T*. For the CSW potential, Q 6 has a non-monotonic behaviour along the 
isotherms with a maximum at pQ max (Fig. EJ). 

The density PQ max for each isotherm lies between pt max and Pt min . In the area between 
PQmax an d Pt min both order parameters decrease for increasing p, hence the liquid gets more 
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FIG. 9: Particle distribution for the first and second peaks of the g(r), normalized by the area of 
each peak, at T* = 0.5: for A = 30 at (a) p* = 0.15 and (b) p* = 0.25, and for A = 500 at (c) 
p* = 0.15 and (d) p* = 0.25. The distance between the maxima of the first and the second peak 
decreases for increasing p and for increasing A. 

disordered with increasing density. This behaviour defines the structural anomaly region 

(PQmax < P < Pt m J- 

By increasing the value of A we find that the region of structural anomaly does not 
contract, but tends assimptotically to a fixed region in the T* — p* plane (Fig. [3]). This weak 
dependence on A suggests that the occurrence of the structural anomaly does not disappear 
for very steep soft-core potentials. This prediction is consistent with the s 2 calculations that 
allow de Oliveira et al.— to argue that the structural anomaly should be observable also for 
the DSW potential, here considered as the limit of the CSW for A — > oo. 



VI. DISCUSSION 



To get more insight into the origin of the anomalies, we study the relation between the 
changes of the coordination shells of the radial distribution function g{r) and the appearence 
of the anomalies. First, we perform a systematic study of the first two peaks of g(r). For 
example, for A = 30 we find that at densities within the LDL phase (Fig. [9^l) the distribution 
of particles of the first peak is maximum at T\ja ~ 1.05, next to the hard-core. For the 
second peak, the distribution is maximum at r 2 /a ~ 1.95, at about the minimum of the 
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atractive well. At densities within the HDL phase many particles overcome the shoulder 
approaching the hard-core. For p* = 0.25 (Fig. [9b ) the maximum of the first peak increases 
without changing its position, at r\/a ~ 1.05, while the distribution of the second peak is 
slightly flattened with a maximum around r/a ~ 1.75. This change by increasing density 
shows that particles in the attractive well of the potential accumulate close to the soft-core 
distance, giving rise to an asymetric distribution of particles around the minimum of the 
attractive well. 

For the highest A considered here, A = 500, we find a similar result at the LDL density 
(Fig. IHt), but at the density p* = 0.25 (Fig. (HH), approaching the LDL-HDL critical point, 
the distribution of particles at the second peak— displays a flat region in the range 1.7 < 
r/a < 2.0 and a maximum around the soft-core distance at Rr/ci = 1.6. This is consistent 
with the fact that for A = 500 the potential has a more penetrable soft-core for r/a > 1.6 
with respect to the case A = 30. As a consequence, the effective soft-core distance for the 
potential with A = 500 is smaller than that for the potential with A = 30. 

This qualitative observation gives rise to a quantitative criterion for the prediction of the 
occurrence of the density and diffusion anomaly when we generalize what has been proposed 
by Yan et al. for potentials with a repulsive ramp^. They observe that, if a and b are 
respectively the starting and ending distance of the repulsive ramp, the anomalies of density 
and diffusion occur only when 0.5 < a/b < 6/7 = 0.86. This is consistent with the fact that 
for real water in standard conditions the ratio between the distances of the first two peaks 
of the oxygen-oxygen g(r) is about 0.6, while in normal liquids it is about 0.5^. 

A straightforward application of the the Yan's criterion to our case is not obvious. One 
could assume that the relevant ratio in the CSW potential is a/R^ = 0.5, or as an alternative 
Rr/Ra = 0.8. With both choices one would conclude that Yan's criterion does not hold for 
the CSW potential because it would always predict the absence (with the first choice) or 
the occurrence (with the second choice) of the anomalies independently on A. 

Hence, to generalize the Yan's criterion to cases in which the effective soft-core distance 
changes with density or A, as in the present case, we adopt as characteristic length-scales 
the maxima of the first two peaks of the g(r) within the HDL phase and define the ratio 
A(A) = r 1 (p)/r2(p)\ PKDL at a temperature below, or at about, the LDL-HDL critical point. 
The ratio A is larger for more penetrable potentials, i. e. for larger A in the case of the 
CSW potential. 
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With this generalized definition, we conclude that the empirical criterion proposed by 
Yan et al.— is valid also for all the cases considered in our work. For the cases in Fig. [9] it 
is, for example, A(A = 30) = 0.55 < A(A = 500) ~ 0.66, where the approximation symbol 
is used because for A = 500 we can evaluate A only at densities approaching Phdl{T).— We 
observe that by increasing A the regions in the plane (T — p) where the density and diffusion 
anomalies occur reduce in a sensible way, while the region of structural anomalies is almost 
unaffected by the change of A (Fig. [3]). 

The considerations above are possibly consistent with the idea that the structure within 
the first two coordination shells is controlling the occurrence of the anomalies. In the 
following we will show that this concept is not leading to any general criterion for the 
occurrence of the anomalies and that the anomalies are, instead, regulated by the structure 
at a larger scale. 

To this goal, we test an idea proposed in Ref.— suggesting that a necessary condition 
for the appearence of the anomalies is a transference of particles from the second shell of 
the fluid to the first shell under isothermal compression. According to this criterion all the 
anomalies appear for the range of p and T where 

dg(r 



t-t _ 9g(r) 

^1,2 : 



dp 
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dp 



< , (9) 
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i. e. Eq.()9]) is a necessary condition. On the other hand, if ITi i2 > 0, no anomalies should 
be observed, i. e. the condition in Eq.Q is also sufficient.— 

To study the implications of Eq. ([9]), we analyse the behaviour of the CSW g{r) upon 
isotermal compression (Fig. [10] and Fig. [TTj) . For all the considered values of A and p* 
we find that it is always n 12 < 0, even for T and p outside the regions of the anomalies. 
Therefore, the condition in Eq. could be necessary, but not sufficient, for the appearance 
of the anomalies. This observation is consistent with the result for the DSW potential, where 
Eq. ([9]) is satisfied (see Fig. 17 in Ref.—) but no density anomaly is observed^"—. 

To understand better why Eq. ([9]) is not a sufficient condition for the anomalies, we study 
how n 12 changes with p and T (Fig. IT2|) . For each T, we find that n 12 (p) is not monotonic, 
consistent with the fact that n 12 (p) is expected to be negative around the anomaly regions 
and to be positive at smaller and larger p. We observe that n 12 (p) reaches a lower minimum 
for lower T, i. e. ITi i2 is more negative where the anomalies are stronger. However, away 
from its minimum ITi j2 (T) is not a monotonic function of T at constant p. Moreover, we find 

17 




FIG. 10: Change of the radial distribution function g(r) for A = 30 with p* and T*. The g(r) at 
densities 0.15 < p* < 0.27 for (a) T* = 0.5 and (b) T* = 0.7. (c) The height of the first peak of 
g(r), i. e. g(ri), always increases for increasing p*, but has a more gradual increase at higher T. 
(d) The height of the second peak of g(r), i. e. gfa), always decreases for increasing p* with a 
weak dependence on T. 

that n 12 (p) is not symmetric with respect to the center of the anomalous regions. Instead, 
IT12 reaches its minimum close to the maximum density with diffusion anomaly (Fig. [T2l) . 
We also observe that for there is no clear threshold below which an anomaly appears, 
contrary to the case of £2 derived from the excess entropy. All these observations suggest 
that there is no condition on n 12 that is sufficient for the occurrence of the anomalies. More 
generally, these results suggest also that the onset of the anomalies is related not only to 
the behavior of the first and second coordination shell of the g(r), but also to the behavior 
of the higher coordination shells, as recently observed by Krekelberg et al. for other models 
of anomalous liquids, including water—. 

To better understand the relevance of the outer shells, following Ref. 6Q we perform an 
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FIG. 11: Radial distribution function g{r) for A = 500. For T* = 0.5 (a) and for T* = 0.7 (b) the 
behaviour of g(r) at different p* is comparable to the case in FigllO[ the main difference being the 
appearance of a fine structure in the first two peaks at high p* , resembling the g(r) of the DSW 



potential 41 . 



analysis of the excess entropy shell by shell, but instead of calculating only the cumulative 
order integral, as in Ref.— , 

rr 

s 2 (r) = -2np / [g(r') In g(r') - g(r') + l]r' 2 dr' (10) 
Jo 

here we consider also its derivative 

dso(r) 



S 2 r 



(11) 

T 



<9 hip 

where S2{r) — > S2 and S 2 (r) — > S 2 for r — > oo. Our analysis shows that S 2 (r) presents many 
oscillations within each coordination shell (Fig. [131 Fig. [HI) . In several cases different shells 
contribute with a different sign to the total S 2 and in some cases the relevant contribution 
to the large-r value of S 2 (r) comes from the third or higher coordinations shell. 
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Density p* 



FIG. 12: Product II12 in Eq.Q as a function of density for A = 30 and different temperatures. 
Symbols mark the lowest and the highest density where at each T a specific anomaly is observed: 
▲ for the maximum Qq and T for the minimum t, marking the structural anomaly region, ■ for 
the limits of the diffusion anomaly region, Q for the limits of the density anomaly region. For other 
values of A, not shown here, we find the same qualitative behavior. 

For example, in the case with A = 30 we observe that at low density (Fig. IT3h.c for 
T* = 0.5 and 0.7, respectively) the major contribution to the structural anomaly comes 
from the first and the second shell, although for p* = 0.17 at T* = 0.7 (Fig. IT3"b ) it is only 
after the forth shell that £2(7") stops to oscillate around the threshold for the structural 
anolmaly. On the other hand, the relevance of the third and higher coordination shells is 
very evident at higher density (Fig. [I3"b.d). where the sign of £2(7") changes in some cases 
within the third shell (p* = 0.24 at T* = 0.5, p* = 0.23 and p* = 24 at T* = 0.7). 

Also for the diffusion and the density anomaly, the relevance of the outer shells is evident, 
especially at high p. For example, in Fig. [T3b there are several cases where the thresholds 
for diffusion and density anomaly are reached only within the third shell. While in the 
majority of cases £2(7") saturates within the first fifth shells, we observe also cases where 
the sixth and possibly higher coordination shells are relevant for the final value of £2, e. g. 
for p* = 0.24 and p* = 0.25 at T* = 0.5 (Fig. fl~3b) with respect to the empirical thresholds 
for the density anomaly and the diffusion anomaly, respectively. Therefore, in general is not 
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FIG. 13: The cumulative order derivative ^(r) calculated for the case A = 30, for different 
temperatures and densities, (a) For T* = 0.5 and density going from p* = 0.16 (lower line) to 
p* = 0.19 (upper line), (b) For T* = 0.5 and density going from p* = 0.20 (upper line) to p* = 0.26 
(lower line), (c) For T* = 0.7 and density going from p* = 0.16 (lower line) to p* = 0.20 (upper 
line), (b) For T* = 0.7 and density going from p* = 0.21 (upper line) to p* = 0.26 (second lower 
line). The values of p* are indicated by the labels near the calculated lines. Horizontal continuous 
lines with labels st, D and p mark the thresholds for structural anomaly, diffusion anomaly, and 
density anomaly, respectively. Vertical dashed lines mark the coordination shells as defined by the 
minima of g(r) in Fig. [TOb.b. 
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possible to establish the occurrence of the anomalies on the base of the local structure (first 
and second shell), but it is necessary the structural information at much larger scale (forth 
and outer shells). 

From Fig. IT3"b.d we observe that the empirical thresholds for diffusion anomaly and density 
anomaly fail to predict correctly the absence of these anomalies at T* = 0.7. The comparison 
with Fig. [3b, indeed, shows that these anolmalies are not found at T* = 0.7, while the 
calculation of £2(7")) over more than two coordination shells, would predict the occurrence 
of the diffusion anomaly for 0.17 < p* < 0.23 and of the density anomaly for 0.18 < p* < 0.22 
(Fig. IT3"b.d). Nevertheless, note that the scale in Fig. [T3"b.d is 40% smaller than in Fig. fT3"a.b 
and £2(7") never exceeds 1.5. 

The same analysis for A = 500 presents even more striking oscillations of £ 2 (r) within 
the first two coordination shells (Fig. [T4"j) . This is, for example, evident for for p* = 0.19 
and T* = 0.7 (Fig. Hlb ) at r/a = 1.5, where the sudden change of slope of the penetrable 
soft-core gives rise to a large negative contribution to £2(7*). However, the contribution of 
the outer shells leads to positive £2. At the same T* and higher p*, similar discontinuities 
at r/a = 1.5 lead to negative £2, as for p* = 0.25 (Fig. [I4"H). A remarkably large positive 
contribution to £ 2 at r/a = 1.5 is observed at T* = 0.5 and p* = 0.23 (Fig. [Tib), in 
correspondence to the limit of the region of structural anomaly (Fig. [3H). However, also in 
this case the contribution coming from the outer shells is comparable to that coming from 
the second shell. 

For A = 500 the system displays structural anomaly at any T, but diffusion and density 
anomalies are present only below T* = 0.5 (Fig. |3H). Nevertheless, from Fig. [THwe observe 
that the empirical thresholds of £2 for diffusion and density anomalies fail at any T, as seen 
for A = 30 (Fig. [T3|) at higher T only. Also for A = 500 the higher the T the lower the 
saturation value of T* 2 {r), as for A = 30. 

Finally, we observe that in all the cases we have analyzed, when £2(7") saturates to a value 
lower than 2.5, both diffusion and density anomalies are not observed. The only possible 
exception is observed for A = 500 at T* = 0.5 and p* = 0.20 (Fig. flla). where £ 2 = 4. 
However, this state point is at the border of the region where both diffusion and density are 
anomalous (Fig. |3H). 
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FIG. 14: As in Fig. H3J but for A = 500. (a) For T* = 0.5 and density going from p* = 0.16 
(lower line) to p* = 0.20 (upper line), (b) For T* = 0.5 and density going from p* = 0.21 (upper 
line) to p* = 0.26 (lower line), (c) For T* = 0.7 and density going from p* = 0.16 (lower line) to 
p* = 0.20 (upper line), (b) For T* = 0.7 and density going from p* = 0.21 (upper line) to p* = 0.25 
(lower line). Vertical dashed lines mark the coordination shells as defined by the minima of g(r) 
in Fig. [IIJt,b. 

VII. SUMMARY AND CONCLUSIONS 

We perform a systematic analysis of the soft-core CSW potential in Eq. (1) by varying 
the slope A of the repulsive shoulder. We consider values from A = 15 (smooth and less 
penetrable at large distance )^ 4 to A = 500 (almost discontinuous and more penetrable at 
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large distance), we calculate their phase diagrams and the occurrence of the anomalies, and 
we compare with the results for the discount inuous (DSW) version of the potential^. 

All the phase diagrams for different values of A display at low density the gas-liquid first 
order phase transition, ending in a critical point C\, and at higher density the LDL-HDL 
first order phase transition, ending in a critical point C<i- By increasing A, both critical 
points aproach the values of the two critical points for the DSW potential (Table I). In this 
sense, the DSW can be considered as a limiting case of the CSW potential. 

For any A we observe the occurrence of the density anomaly, diffusion anomaly, and 
structural anomaly with the same hierarchy as for water-^, i. e. the density anomaly 
occurrs in the inner region and the structural anomaly in the most external region. The 
difusion anomaly region always includes the density anomaly region, being the two regions 
only slightly separated, with decreasing separation for decreasing p and increasing A (Fig. [3]). 

The region of structural anomaly is always present and almost independent of the value 
of A. It slightly contracts toward a limiting finite region for increasing A. On the other 
hand, the regions of diffusion and density anomalies narrow for increasing A. A similar 
behavior has been found recently for purely repulsive soft-core potentials 61 . Our finding 
gives a rationale for the absence of density anomaly for the DSW potential^"—, here seen 
as the limit of the CSW for A — > oo. 

We find that the excess entropy analysis gives an approximate estimate of the anomalous 
regions. However, comparison with Ref.— shows that the behaviour of £2 for the DSW 
potential is still far from what we obtain for A = 500, consistent with a slow convergence 
of the CSW toward the DSW potential when A is increased. 

We tested some ideas recently suggested as possible explanations for the anomalous be- 
havior of other soft-core potentials. First, we generalized to our CSW potential the Yan 
et al. 58 criterion for the occurence of the diffusion and density anomalies in ramp-like po- 
tentials with two characteristic length-scales. To this goal we propose to use as effective 
length-scales the distances of maxima of the first and second peak of g(r) within the HDL 
phase. While Yan et al. criterion does not apply to our CSW potentials, our generalized 
criterion is consistent with our data for different values of A and is possibly valid for other 
continuous potentials with competing length-scales. 

Next, we observe that the condition in Eq. (Q proposed by de Oliveira et al.— is not 
enough to mark the occurence of the anomalies. This condition relies on the structural 
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changes of the first and second coordination shell, while from our calculations we conclude 
that the occurrence of anomalies is related to long range structural changes. To show this 
result we consider the derivative S2( r ), Eq. (jlip . of the cumulative order integral introduced 
by Krekelberg et al. 60 . The study of S 2 (r) shows that the anomalous behavior of the 
potential is dominated at high p and low T by the structural changes up to the fourth 
coordination shell, and possibly by even outer shells. We observe that by increasing A 
(more penetrable and sharp soft-core) the contribution to S2 from the first and second shell 
becomes more relevant. We conclude that the gradual disappearing of the anomalies for 
large A is regulated by the appearence of a more structured short-range order (first and 
second coordination shell) that dominates over the effect of the larger-range structure. 
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